Methods of Energy Transport 


There are 3 methods for energy transport within a star 


e Radiative absorption and emission 
e Convection 
e Conduction (collisions between particles) 


Conduction is only important under degenerate conditions (i.e., 
in white dwarfs and neutron stars). However, in normal stars, 
both radiation and convection can be important. 


Radiative Transport of Energy 


Consider a system of particles diffusing across a boundary in the 
z direction. Let 


n = the particle density 

v = the mean velocity of the particles 

l = the particle mean free path. For particles, l is defined as 
the mean distance a particle will move before its cross- 
section for collision, gø, encounters another particle, ol = 
1/n 


If we assume that the motions of the particles are isotropic, then 
on average, ~ 1/3 of the particles will have their motions primar- 
ily in the z direction, and ~ 1/2 of those will be moving toward 
the boundary, as opposed to away from it. Thus, the flux of 
particles diffusing across the boundary will be 


1 
Fi = g Al (3.1.1) 


Similarly, the flux of particles diffusing across the boundary from 
the other direction is 


1 
Fy = g 7H Uz+l (3.1.2) 


If vz—1 & vz41, then the net flux of particles will be 
1 
fi — Fo = ae [n _ Nz+1| 


If the mean free path is smaller than the density gradient, this is 
simply 


1 
F= ——y | — F=-—D ae 
ga — Vn (3:13) 


where the diffusion coefficient, D = 40l. (This is Fick’s law.) 


To compute the flux of radiative energy across a boundary, treat 
the photons as particles, and recall that the energy density of 
radiation, U, is 


U=aT* 
ei 1 dU 1 dT 
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3° Ph Gp 5 eha (4a =) 
For photons, the mean free path between collisions is 
1 
bun — 
ph Kp 


where « is the mean absorption coefficient (with units of area per 
unit mass), so the diffusion equation can be written 


AacT? dT 
3kp dr 


(3.1.4) 


To put this in terms of an expression useful for stellar structure 
calculations, first substitute £ for F, 


L AacT? dT 


Anr2 — 3kp dr 
and then change the equation to the Lagrangian form 


dT 3KL 


ge ee ally 
dM 647r2acr4T? (3.1.5) 


We can now put this in terms of the radiative temperature gra- 
dient for a star in (at least approximate) hydrostatic equilibrium 


YV dla T = Madar /M dP 
ed \dnPjJ a TdM/ PdM 


Substituting the Lagrangian equation for the radiative tempera- 
ture gradient (3.1.5) and the pressure (2.2.4) we get 


3KLP 


Vrad = Torac MTI 


(3.1.6) 


and the radiative plus conductive flux at any point in the star is 


4acGMT* 
3K 72 P 


Poe Va (3.1.7) 


Note that the diffusion approximation is only good when the 
photon mean free path is short compared to the scale of the 
density gradient. Near the surface, the approximation breaks 
down. 


Criteria for Convection 


To decide whether convection occurs in a star, begin by consid- 
ering a blob of material with volume V in pressure equilibrium 
with its surroundings. Initially, both the blob and its surround- 
ings have temperature T) and density pı. Now suppose the blob’s 
temperature is perturbed, such that its new temperature is T}. 
If Ti > T, and pressure balance is maintained, then the equation 
of state of the gas demands that p’ < p. Given a local gravity, g, 
the blob will feel a buoyance force per unit volume 


pig — pig 
As a result, the blob will rise a distance Ar, to a location in the 
star with Tə and p2. As it does so, the blob will expand to a larger 
volume (to maintain pressure equilibrium with its surroundings) 
and its new density will be p5. The buoyancy force on the blob 
will now be 

p29 — P29 
If p < p2, then the blob will continue to rise, and the star will 


be dynamically unstable to convection; if pS > p2 then the blob 
will sink back to its original position. In other words, convection 


will not occur if p 7 
p p 

— — 3.2.1 
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where the subscript i denotes how the internal density of the 
blob has changed with radius, and the subscript s indicates the 
gradient in the star. 


To translate this condition for stability into a more tractable 
form, begin by using the equation of state, p = p(P,T, u). dp is 
then related to P, T, and u by 


_ (2 a) ap (8 
dp = (55 ) ap (2) ar+ (22) du 


SO 
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where a is the isothermal compressibility, 6 is the volume coef- 
ficient of expansion, and ọ is the chemical potential coefficient. 
Thus, the condition for stability is then 


adP dT Y du adP dT Y du 
Pd Tdr udr Par Tdr udr 


Since, by definition, the blob must remain in pressure equilib- 
rium with its surroundings, the first terms on either side of the 
inequality cancel. Furthermore, since the blob does not change 
its chemical composition while moving, (du/dr); = 0. This leaves 


ô dT > ô dr pdu 

Tdr); Tdr pudrj, 
Now if we multiply both sides of the inequality by the pressure 
scale height, 


dP dr 
Ap=-—P We ae (3.2.3) 


the condition for stability becomes 
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Note that if the star is stable, the term on the left is Vraq, the 
star’s radiative temperature gradient. Moreover, the first term 
on the right is nothing more than the adiabatic temperature gra- 


dient. Thus, for a chemically homogeneous region of a star, the 
condition for stability is simply 


Vead < Vad (3.2.5) 


This is the Schwartzschild criterion for dynamical stability. If 
the star is chemically inhomogeneous, then the Ledoux stability 
criterion must be used: 


p { dinu 
ae 
Vrad < Vaa a 5 (5-5) (3 6) 


Convective Transport of Energy 


In a radiative region of a star, the temperature gradient is given 
by Vraa. However, energy transport by radiation (and conduc- 
tion) is not always efficient. In that case, convection transport of 
energy will occur. By definition, if radiation and conduction are 
not working, then energy is not being lost from the individual 
convective cells, and the stellar temperature gradient V = Vag. 
If energy transport by radiation and convection are comparable 
then Vraa > V > Vad. 


In practice, when convection occurs in stellar cores, V ~ Vaq. 
Only in the outer parts of some stars do superadiabatic conditions 
exist. 


